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Abstract: Six semi-detached Algol-type binaries lacking a period analysis were chosen to test 
for a presence of a third body. The O — C diagrams of these binaries were analyzed with 
the least-squares method by using all available times of minima. Also fourteen new minima, 
obtained from our observations, were included in the present research. The light-time effect was 
adopted as a main factor for the detailed description of the long-term period changes. Third 
bodies were found with orbital periods from 46 up to 84 years, and eccentricities from 0.0 to 
0.78 for the selected binaries. The mass functions and the minimal masses of such bodies were 
also calculated. 

Keywords: stars: binaries: eclipsing; stars: individual: DK Cep, TY Del, RR Dra, TZ Eri, 
VX Lac, UZ Sge; stars: fundamental parameters 

1 Introduction 

Eclipsing Binaries (hereafter EBs) are excellent objects for determining the physical properties 
of stars and detecting additional components in these systems. The long-time behavior of the 
period of an EE could reveal the presence of another component orbiting with the EB around 
the common center of mass. Photometric observations of EBs sometimes cover more than a 
century, so the periods of such third bodies could reach the same time interval. 

The motion around the barycenter causes apparent changes of the observed binary period 
with a period corresponding to the orbital pe riod o f the t hird body, called the Light-Time 
Effect (or ' li ght-tr avel time', hereafter LITE). Ilrwini (119591 ) improved the method developed 



by IWoltjen (119221 ) for analyzing the long-term variation of the times of minima caused by 
the third body orbiting the e c lipsing pair. Useful c o mmen ts and limitations were discussed 
by Frieboes-Conde fc Herczeg ( 1973 ) and by Mayer ( 1990l ). Nowadays there are more than 
one hundred EBs showing L I TE, where the e f fect i s certainly p r esented or supposed ( s ee e.g . 
Borkovits fc HegeduesI (ll99fih . lAlbavrak et all (|l999h . IWolf et (l2004h . lHoffman et al.l (l2006h . 
etc.). For the apparent peri o d cha nges in EBs according to their O — C diagrams, see e.g. the 
catalogue by iKreiner et al.l (120011 ). From the same paper the look of the O — C diagrams, 
presented in the present paper, was adopted. In the figures [T] to [H] of the present work the full 
circles represent the primaries and the open circles the secondaries, the bigger the point, the 
bigger t h e weig ht. For the limitations and consequences of the O — C diagram analysis, see e.g. 
Sterkenl tOQ^ . 



^Email: zasche@sirrah.troja.mfF.cuni.cz 
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The computation of the parameters of the third-body orbit is a classical inverse problem 
with 5 parameters to be found - p^, Tq, A, u, 63, which indicate the period of the third body, the 
periastron passage, the semi-amplitude of the light-time effect, the argumen t of periastron and 
the eccentricity, respectively (for a detailed description see e.g. Mayer 1990l ). The ephemerides 
for the individual system [JDq and P for the linear one and q for the quadratic one) have to be 
calculated together with the parameters of LITE. The mass function /(M3) and the minimal 
mass of the third component M^^min = ■ sin ^3 (for = 90°) could be computed from this set 
of parameters. The weights assigned to individual observations were used as following: w = 1 
for visual observations, 3-5 for photographic and 10 for CCD and photoelectric observations. 
The computing code itself could be downloaded via the web pages of the authoiEI. 

All of the selected systems are Algol-type EBs, semi-detached and also their spectral types 
are similar. The primary components have spectral types from B to F, while the secondaries 
from G to K. Except for 14 observations in Table [TJ all the times of minima used in this 
paper were collected from the published literature and from minima databases available in the 
internet. 

According to a recent paper on the period changes in Algols by lHoffman et al.l (l2006l ). there 
could be a connection between the spectral type of the secondary component and the nature 
of the period changes. Systems with spectral types of secondaries later than F5 show O — C 
variations, which could be c aused by t h e magnetic activi t y cycles and convec tive envelopes. 
This effect was discussed by iHalll (jl989l ). lApplegatd (119921 ). iLanza et al.l (119981 ). etc. The role 
of magnetic cycles on the period changes is discussed below, but due to lack of information 
about the systems such an analysis is a difficult task. For some of the systems selected in 
this paper the spectral types of secondaries are known with a low confidence level, light-curve 
analysis is missing and spectroscopy has never been done. 



2 Observations 

The photometric observations were obtained by CCD detectors at two different observatories. 

The first one is situated in Ondfejov Observatory, Czech Republic. The 65-cm telescope 
was equipped with the Apogee AP-7 CCD camera. All of the observations were secured during 
2006 and 2007 and only R filter was used. The exposure times depend on the particular target 
and observing conditions, ranging from 20 to 60 seconds. 

The second observatory is situated at Athens University campus, Athens, Greece. The 40- 
cm telescope is equipped with the SBIG ST-8XMEI CCD camera. All of the observations were 
obtained during 2007, using only the R filter. The exposure times depend on the individual 
stars and observing conditions, ranging from 15 to 120 seconds. 

^ http: / / sirrah.troja.mff.cuni.cz/~zasche/ 
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3 Analysis of individual systems 



3.1 DK Cep 



The first system is the eclipsing binary DK Cep. This neglected system is about 12.2 mag 
bright in B filter. The spectrum was classified as A8 + G4IV in ISvechnikov &: Kuznetsova 



(1l990l ). but only on the basis of pho tometric indices. Its photometric variability was discovered 
in 1966 by iRomano fc Perissinottd . but since then no analysis of this system was performed 
(neither photometric nor spectroscopic one). Also period analysis has not been performed so 
far. 

Altogether 157 observations of minimum light were collected for the analysis. One secondary 
and two primary times of minima were also observed at Athens observatory (see Table [I]). 

The analysis was carried out and the resultant plot is given in FiglH where the theoretical 
curve is shown together with the individual observations. The parameters computed from our 
analysis are given in Table [21 Regrettably, the phase near the periastron is not covered with 
the data. This affects the values of u and e and their respective errors, so the orbit could 
have slightly different eccentricity. This hypothesis could be proved in the upcoming periastron 
passage, which will occur near 2010. 

From the parameters of the third-body orbit, the predicted minimal mass of this component 
can be also derived. Using the masses of the primary and secondary of th e EB pair as Mi = 
1.65 Mq and M2 = 0.92 Mq (according to ISvechnikov fc Kuznetsovalll990h . the minimal mass 
results in Ms^rnin = 0.32 M0. Due to large difference between this value and the values of 
masses of the eclipsing pair components, one could speculate about the pos sible detectabilit y 
of such a body. This mass indicates a spectral type about M3 (according to iHarmaned Il988l ) . 



Table 1: New times of minima based on CCD observations (Kwee-van Woerden (119561 ) method 



was used). All of them were obtained in R filer, 
abbreviations see the text. 



For the explanation of the observatory 
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Figure 1: The O — C diagram of DK Cep. The individual observations are shown as dots (primary) and open 
circles (secondary), the small ones for visual and the large ones for CCD and photoelectric observations, the 
bigger the symbol, the bigger the weight. The curve represents the predicted LITE variation (see the text). 



Table 2: The final results (part 1): Parameters of the third-body orbits from the analysis of 
times of minima. The table is divided into two parts. In the first one the computed parameters 
are presented, while in the second one the derived quantities are given (the mass M12 = M1 + M2 
is taken from the literature). 



Parameter 


DK Cep 


TY Del 


RR Dra 


JDq [HJD] 


2433590.556 ± 0.005 


2442959.471 ± 0.001 


2434913.728 ± 0.022 


P [day] 


0.9859085 ± 0.0000004 


1.1911264 ±0.0000002 


2.8312140 ± 0.0000053 


P3 [yr] 


31.3 ±2.5 


64.9 ±2.3 


84.3 ±0.6 


To [HJD] 


2454900 ± 1200 


2449192 ± 1017 


2450058 ± 442 


uj [deg] 


273 ± 32 


38.5 ±18.2 


110.2 ±3.5 


e 


0.780 ± 0.256 


0.221 ±0.055 


0.503 ± 0.028 


A [day] 


0.009 ± 0.002 


0.0268 ± 0.0012 


0.0731 ± 0.0017 


q [day] 






(-126.2 ±0.1) • 10-1" 


M12 [ Mq] 


2.57 


3.64 


2.75 


/(M3) [ Mq] 


0.0039 ± 0.0002 


0.025 ± 0.001 


0.300 ± 0.002 


M3,„i„ [ Mq] 


0.32 ±0.01 


0.79 ± 0.07 


1.85 ±0.09 



The magnitude difference between the primary and the tertiary component is therefore more 
than 6 mag, or the amount of the third light in the light curve less than 1 percent, which 
is undetectable. Also, in the spectrum of the system such a body will be, probably, hardly 
detectable. Since the distance to the system is not known, one cannot estimate the predicted 
angular separation of the third component. 



3.2 TY Del 

TY Del (AN 141.1935, BD+12 4539) is another EB with apparently variable period. The 
spectrum was classified as B9+G0IV (IHoffman et al.l 120061 ) and the relative brightness about 
10.1 mag in V filter. There is a consensus about t he spectral types of th e comp onents of 
TY Del, but there is a difference between the masses. iBrancewicz fc DworakI (jl980h . and also 
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Figure 2: The O — C diagram of TY Del (for the description see FigdJ. 



Buddind fll984Jl and iBudding et al.l (120041) have proposed masses Mi = 5 Mq, M2 = 2 M 



while ISvechnikov &: Kuznetsoval (Il990l ) gave Mi 
hkely to the proposed spectral types. 



0) 

2.8 Mq, M2 = 0.84 Mf,,. The latter are more 



0- 



The star was discovered to be a variable by iHoffmeisterl (119351). There was on ly one attempt 
to observe the whole light curve of TY Del photoelectrically by iFaulknerl (119831 ). unfortunately 
only about a half of the curve w as obs e rved. No analysis of these data was carried out so 
far. The star was also studied by ICookl (119931 ) on the basis of his visual observations for the 
long-time scale intrinsic variations, but the results were not very conclusive. 

The spectroscopic observations in Ha were obtained by I Vesper et al.l (I2OOII ). They concluded 
that there was no activity in Ha and no evidence for the mass transfer structures was found in 
this system. 

Altogether 368 times of minima were collected, from which only 5 were omitted due to their 
large scatter. One period of the third body is already sufficiently covered by the data points (see 
Fig. [2]), but further observations are still needed. From the LITE para meters (see Table [2)) and 
with the approximate masses of the EB components as M12 = 3.64 (ISvechnikov fc Kuznetsova 



19901 ) it was possible to derive the minimal mass of the third component as M^^min = 0.67 Mq. 
Due to lack of any other observations, this hypothesis cannot be proved. The spectral types 
and masses were derived only on the basis of the photometric indices and are not very con- 
clusive. A spectroscopic analysis, as well as an analysis of the light curve of the system is 
needed, but the third light is hardly detectable in the light-curve solution. Regrettably, the 
star was not measured by Hipparcos, so the distance is not known and one cannot derive the 
predicted angular separation of the third component. There is also clearly visible some ad- 
ditional short-periodic variation, which could not be described by the LITE. For a discussion 
about its possible explanation see Section 13.71 



3.3 RR Dra 

Another eclipsing binary which exhibits apparent period changes is RR Dra (AN 188.1904). Its 
relative brightness is about 9.8 mag in V filter, and the spectrum was classified as A2+G8IV 
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Figure 3: The O — C diagram of RR Dra (for the descriptfon see Fig[T]). 



(ISvechnikoy fc Kuznetsoval. 119901). wh i le lYoon et al.l (119941 ) proposed a later spectral type KO of 
secondary. ISvechnikov &: Kuznetsoval (119901 ) gave the masses M\ = 2.15 M© and M2 = 0.6 M©. 
The primary minimum is very deep, about 3.5 mag and the or bital period is about 2.8 days. 

The star was discovered to be a variable by ICeraskil (1l905l ). The minimum is so deep that 
visual observers could also provide reliable observations. That is the re ason why most of the 
collected times of minima are visual ones (193 out of 219). iKreinerl (jl97ll ) collected all available 
minima for the period analysis. The long-term increase of the period, evident from his O — C 
diagram, is very probably due to the mass tra nsfer between the components. The most recent 
period study of this system was performed by iQian et al.l (120021 ). who considered (besides the 
mass transfer) the abrupt period jumps - altogether 8 jumps were introduced to describe the 
O — C diagram in detail. Almost the same goodness of fit could be reached by applying the 
LITE hypothesis besides the mass transfer. 

A total of 219 times of minima were used for the present analysis. The O — C plot is 
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Figure 5: The O — C diagram of TZ Eri (for the description see Figd]). 



presented in FigJ^l where LITE and the quadratic term were plotted together. In FigjH only 
LITE is shown. It is obvious that the period increase is very rapid, and the amplitude of 
LITE is still quite high. This leads to the relatively high mass function, which results in high 
predicted minimal mass Ms^min = 1-85 Mq. This is larger than the secondary and should be 
probably observable in the light-curve solution as well as in the spectrum. Regrettably, no such 
analysis has been performed so far. 

The quadratic term coefficient q = (126.1 ± 0.1) ■ 10~^° day leads to a period change about 
3.25 ■ 10~^ day/yr. From this value the conservative mass transfer rate could be derived M = 
3.2 ■ 10^"^ M0/yr. This relatively high value of mass transfer rate arises from the very rapid 
period change, which was attributed t o the quadratic epheme ris. The spectroscopic observations 
during the primary eclipse made by iKaitchuck et al.l (Il985l ) indicate a possible presence of a 
transient accretion disc in the system. The presence of such a disc also supports the hypothesis 
of mass transfer in the system. 

For the estimation of the angular separation of the third body and the astrometric confir- 
mation of the LITE hypothesis, the distance to the system has to be known. The star was not 
measured by Hipparcos a nd the distance is n ot known precisely. The only information about 
the distance comes from Kharchenkol ( 2001 ). where is surprisingly introduced an inaccurate 
value of the parallax vr = (0.40 ± 11.50) mas. A distance with such a large error is useless for 
the estimation of the angular separation of the predicted component. 



3.4 TZ Eri 

The system TZ Eri (AN 40.1929, BD-06 880) is an EB with orbital period of about 2.6 days 
and apparent brightness of about 9.7 mag in V filter. It has a deep primary minimum (about 
2.8 mag), so the visual observations could be also reliab le. 

Its variability was discovered by iHoffmeister (jl929h^ who recog nized the system to be of 
Algol- type. The spe ctral type was f irstly classified by ICannonl ( fl93l as F. Later, the spectrum 
was re-classified by iBarblan et al.l (119981 ) as A5/6 V (primary) and KO/1 III (secondary). In 
this latter paper the light-curve observations in the Geneva 7-color photometric system were 
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Table 3: The final results, part 2. The distance d of the system TZ Eri is taken from the 
literature. 
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analyzed together with the radial-velocity curves of both components. The Wilson-Devinney 
code was used, resulting in a set of parameters descri bing both componen ts. For our analysis 



the masses Mi = 1.97 Mq, M2 = 0.37 M© were used f lBarblan et al.l . \l9m 



Several studies about the pres e nce of an accretion disc i n the system h a ye als o been made 



by iKaitchuck fc HonevcuttI fll982f ). iKaitchuk fc ParkI fllQSSf ). IVesper et al.l fl200lf ) and others 



This disc as well as mass transfer from the secondary to the primary is in agreement with our 
result about the increasing orbital period (see be low). The system y yas also included in the 
sample of Algol-type binaries with radio emission f Umana et al.l . 1998 ). Recently, the star was 
investigated acc;ording to a possible connection between the orbital and pulsational periods (see 



Soydugan et al.ll2006l ). 



The analysis of the long-term period changes was based on a set of 108 observations (mostly 
the visual ones). The resultant O — C diagram is presented in FigJ^] and the parameters of 
the predicted LITE are given in Table [3l The minimal mass of the third component results in 
M^,min = 1-3 Mq, ot the spectral type F6 (according to lHarmanedll988l ). Such a body could 
be detected in the light-curve solution as well as in the spectra of TZ Eri. Regrettabl y , ther e 



was no attempt to detect such a body during the detailed analysis by iBarblan et al.l (119981 ) 



The long-term period increase is due to the mass transfer from the secondary to the primary, 
with the conservative mass-transfer rate M = 3.2 ■ 10~^ Mq/jt. 

Despite the fact that the star was not observed by Hipparcos, iBarblan et al.l ( 119981 ) esti- 
mated the photometric distance to d = (270 ± 12) pc. Assuming a coplanar orbit of the third 
component, then M3 = M3 and one could calculate the predicted angular separation of the 
third body to a = 77 mas and a magnitude difference between the eclipsing pair and the pre- 
dicted component about 1.7 mag. Therefore, the final result is that such a body is detectable 
with the modern stellar interferometers. 



3.5 VX Lac 

The eclipsing binary system VX Lac (BD -|-37 4662, GSC 03214-01295) has its apparent bright- 
ness about 10.55 mag in V filter. The orbital period is close to one day. Its spectral type was 
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Figure 6: The O — C diagram of VX Lac (for the description see Fig[T]). 



firstly derived by ICannod ( 
(ISvechnikov fc Kuznetsova 



1990) 



1934 ) as FO, nowadays the adopted spectral type is FO + K4IV 



Unfortuna tely, neither light curve analysis nor spectroscopic mass ratio exists for the system. 
Kreinen ( 1l97ll ) included this star to his study of systems with period changes, but at that time 
the LITE variation could not be identified conclusively. Nowadays there are 370 times of 
minima. It is clearly visible from this data set that, besides the quadratic term, there is a 
periodic variation in the times of minima (see FiglH]). The parameters of the fit are given in 
Table El 

The mass transfer coefficient leads to the conservative mass transfer rat e about M = 4.7 • 
10"^ M0/yr. Such a result is in agreement with the previous study by I Chid (1l998l ). who 
concluded that the star is currently in the mass loss state. 

The LITE parameters indicates the predicted minimal third-body mass about 0.4 Mq. The 
secondary component (M2 = 0.47 Mq) is only slightly more massive, but, despite this fact, 
the expected magnitude difference between the primary and the tertiary component is too high 
(about 5 mag) to be detected. The derived contribution of the third star to the total luminosity 
is only about 1 percent and therefore is hardly detectable. Having no information about its 
distance one cannot estimate the angular separation of the third body. The additional variation 
in O — C diagram after subtracting the LITE and quadratic term is shown in FiglH] (see the 
possible explanation in Section [3Tll . 



3.6 UZ Sge 

The EB system UZ Sge (A N 435.1936) has orbital pe r iod o f about 2.2 days and its spectrum 
is classified as A3V+G0IV (ISvechnikov fc Kuznetsoval . Il990l). 

The photometric variabihty of UZ Sge was discovered by lGuthnick fc Schnelleii (119391 ). Since 
then, there was no attempt for any detailed analysi s, neither photo metric nor spectroscopic 
one. The only spectroscopic observation was done by iHalbedell (119841 ) for the derivation of the 
spectral type of primary component. 

Five new times of minima were obtained for this paper (see Table [T]). Altogether 125 
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Figure 7: The O — C diagram of UZ Sgc (for the description see Figd]). 



measurements of times of minima were found in literature, but 14 observations were neglected 
due to their large scatter . Tak ing as masses of the individual components those given by 
Svechnikov fc Kuznetsova Ml = 2.05 M0 and M2 = 0.29 M©, the minimal mass of the 

third component results in M^^rnin = 0.65 M©. But due to absence of any detailed analysis of 
this system, this value cannot be proved. 



3.7 Alternative explanation of the orbital period changes 

One can see some additional non-periodic variations in the O — C diagrams, which cannot be 
described by applying only the LITE hypothesis. In Fig. M four cases with the most evident 
variations are shown. The amplitudes of these variations are usually about 10 minutes in O — C 
diagram and these are not strictly periodic. This could be caused by the presenc e of stellar 



(1992 


), 


Lanza et al. 


(1998). or 



spectral types of the secondary components are later than F5 (see IZavala et al.l (120021 ) for a 
detailed analysis). To conclude, for a better description of the observed period variations of 
these systems, the magnetic activity cycles could be present together with the LITE. On the 
other hand, one has to take into consideration that the spect ral types of most of these b i naries 
were derived only on the basis of their photometric indices ( Svechnikov Sz Kuznetsova . 1990l ) 
and are not very reliable. 

The only sy stem, where one could estimate the variation of the quadruple moment (see 
Applegatelll992l ) required to explain the long-term period variations, is TZ Eri, where the semi- 
major axis of the orbit from the light curve and radial velocity curve solution is known. Using 

the following equation 

AP = A^/2{1 - cos{27rP/j93}) 

( Rovithis-Livaniou et al.l . 2000l ) one could compute the amphtude of the period oscillation. The 
period variation AP/P can be used for calculating the variation of the quadruple moment AQ, 
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Figure 8: The O — C diagrams after subtraction of LITE. 



using the equation (ILanza fc Rodond |2002| ) 

AP 



P 



-9 



This quantity results in AQ = 3.6- 10^^ g-cm^, which is well within the limits for active binaries 
and therefore the variation in TZ Eri could be also explained by this mechanism. 



4 Discussion and conclusions 

Six Algol-type semi-detached eclipsing binaries were analyzed for the presence of LITE on 
the basis of their O — C diagram analysis and the times-of-minima variations. A few new 
observations of these systems were obtained and used in the present analysis. All of the studied 
systems show apparent changes of their orbital periods, which could be explained as a result of 
orbiting the EB with a third component around their common center of mass. 

Such a variation has usually a period of the order of decades, as one can see from FigsilHTl 
which can be described by applying the LITE hypothesis sufficiently. In three cases (RR Dra, 
TZ Eri and VX Lac) the quadratic term in the light elements was also used. This could be 
explained as a mass transfer between the two components, which is a common procedure in 
semi-detached systems. The conservative mass-transfer rates were calculated. 

About half of the systems show also an additional variation in O — C diagram after the 
subtraction of LITE. Such a variation is not strictly periodic, it has usually an amplitude about 
10 minutes and a "period" from 5 to 20 years. Because all of the systems have the spectral types 
of secondaries later than F5, the variation could be caused by the convection and the magnetic 
activity in these binaries. This explanation would clarify the non-periodicity and the changes 
in amplitude of such variations, as well the existence or non existence of such phenomena in 
some EBs. 

Regrettably, in most of the systems no detailed analysis (neither photometric nor spectro- 
scopic) has been made so far. The spectral types and the masses of the individual components 
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in the systems are only approximate, so the parameters of the predicted third bodies are also 
affected by relatively large errors. Due to missing information about the distances to these 
binaries, any prediction about the angular separation could be done. It is obvious that only 
further detailed photometric, as well as spectroscopic and interferometric analysis would reveal 
the nature of these systems and confirm or reject the third-body hypothesis. 
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